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Abstract 

We apply the Euler-Maclaurin formula to find the asymptotic expansion of 
the sums Y2=i ( lo S k) p /k\ J> 9 (logfc) p , £(log k) p /(n- k) q , £ l/& 9 (log k) p in 
closed form to arbitrary order (p, q G N). The expressions often simplify consid- 
erably and the coefficients are recognizable constants. The constant terms of the 
asymptotics are either (( p \±q) (first two sums), (third sum) or yield novel 
mathematical constants (fourth sum). This allows numerical computation of 
C^(±g) faster than any current software. One of the constants also appears in 
the expansion of the function Y2 n >2( n log^) _s around the singularity at s = 1; 
this requires the asymptotics of the incomplete gamma function. The manipu- 
lations involve polylogs for which we find a representation in terms of Nielsen 
integrals, as well as mysterious conjectures for Bernoulli numbers. Applications 
include the determination of the asymptotic growth of the Taylor coefficients of 
(— zj log(l — z)) k . We also give the asymptotics of Stirling numbers of first kind 
and their formula in terms of harmonic numbers. 

To appear in: Results in Mathematics. 



1 Introduction 



This paper is about concrete mathematics. It gathers several results about asymptotic theory, half of 
which are obtained from the Euler-Maclaurin formula. A few by-products offer themselves, such as the 
asymptotics of the incomplete gamma function, or the study of the complex function ^2 n ^ 2 ( n ^°& n ) s 
with a singularity at s = 1, or representations of polylogs in terms of Nielsen integrals, or properties of 
Stirling numbers, or some identities about Bernoulli numbers. We also summarise three ways of obtaining 
the asymptotic growth of the Taylor coefficients of (— 2/log(l — z)) k . Much of the contents may not be 
new - let alone ground-breaking, but the interest of the paper lies in the way all these objects tie the knot 
and pop up by studying a few simple problems; it will offer some surprises to the curious and hands-on 
mathematician. 

To begin with, we recall the Euler-Maclaurin formula: 

n-1 „„ m „ 

£/(*)= / f{x)dx-l\f{n)- m + Y.T^lf^-'Hn)- f^ k -'\a)]+errov 

k=a Ja fc=l [ZK >- 

where the error term is 0( ( 27r Um ) J" \f^ 2m \x)\dx. The values of are j^, — j^q, 30I40' 

We shall be interested in the limit of large n, keeping a fixed. When ordering the terms in decreasing 
orders of n, the quantity — Y^k=i (fty!^ 2 *^^ ( a ) wm contribute to the constant term. The constant term 
will be exact when all orders have been taken into account (ie. m — ► 00). Since this means adding always 
bigger chunks (£>2fc ~ (2k) 2k ), we would end up with an infinite value for the constant term. In practice, 
the exact value of the constant term has to be computed from another approach. However, the formal 

infinite sum involving Bernoulli numbers appears most useful, as it behaves linearly: adding two such 

(2) 

sums (from the asymptotics of H n and iJ„ ' , say) will stand for a constant whose exact value is the sum 
of the two exact values of the respective constant terms. 

We shall use this trick in section 2 to write down the exact constants hiding behind formal sums. 
They will prove useful in subsequent sections to derive the coefficients in the asymptotic expansions (for 
large n) of the four sums that we consider in sections 4,5,6,7 respectively: 

y^Oog^ y^ (1 , )P y^M y 1 

for p, q e N. We shall write their asymptotics in closed form to arbitrary order of n. In particular, we 
can write down (^(±q) and the Stieltjes constants 7 P as formal sums over rational numbers. In this 
formal sense, 7p = {-1^(^(1). 

The coefficients in the asymptotic expansions often contain Stirling numbers of the first kind, or their 
close relative which we denote by S ryS ,t := ^ (sum over all integers ij such that s ^ i\ < ... < i r ^ t). 

Section 3 expresses these numbers in terms of harmonic numbers, which allows a rapid deduction of their 
asymptotics to arbitrary order. The formula can be inverted to express harmonic numbers in terms of 
Stirling numbers. 
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The asymptotic expansion of the four sums, presented in sections 4,5,6,7, can be easily derived from 
the Euler-Maclaurin formula for the first two sums but involves intricate algebra for the latter two. In 
those cases, the expansion was first found empirically using the asympfc trick (appendix). The coefficients 
are all rational numbers except for the constant terms (C^ (i?) fo r the first two sums, unknown constants 
for the fourth sum). For the third sum, the constant term vanishes but C^ P H~ <?) occurs at higher orders 
(irrational). The asymp^ trick gives us sufficient digits of a coefficient c; we then can use the PARI 
software to find a vanishing integer linear combination of 1, c, C'(7), say, if one suspected there was a £'(7) 
hiding behind c. The proper linear combination requires often guesswork. As an application, knowing 
a large number of terms of the asymptotic expansion of the first two sums allows one to compute the 
constant term (±g) to arbitrary precision more rapidly than any current mathematical software; the 
asympfc trick can also enhance speed. 

As an application, we derive in section 8 the asymptotic growth of the coefficients in the Taylor 
expansion of (— z/log(l — z)) , via a convolution from the ansatz at k = 1 (the latter known to Polya). 
The result appeared in two other contexts ( [N-61] and [FO-90]) which we recapitulate for the interested 
reader. 

As advertised, the fourth sum (section 7) gives birth to a 2d- array of unknown mathematical constants, 
C PA , that converge to the values of l/(2 9 (log2) p ) when p,q — ► oo; only C^o and Ci t i have appeared 
(indirectly) before in the literature. Section 9 verifies that C^i, which occurs in X)fc=i k log k ~ logl°g n + 
Ci.i + 0( - il gn ), also occurs in the constant term of the asymptotic expansion of the following complex 
function around its singularity at s = 1: 

E (nl0 g n) , * - lQ g(* - 1) + - 7 + 0{(s I) log 2 (s - 1)) as 

This involves the asymptotics of the incomplete gamma function. 

In order to prove the asymptotic expansion of the third sum (section 6) via the Euler-Maclaurin 
formula, one needs to track down surprising cancellations. The manipulations involve a particular repre- 
sentation of polylogs by Nielsen integrals &i. p (x) presented in section 10 : 

Uj(l - x) = £ (C(j - r) - S^^fr)) ^ 1 - 30 , 

r=0 

wherein the term with £(1) should be dropped. The generalised polylogs Li slv .. iSji (a;) := *i 1 a k (sum 

n 1 ...n fc 

over integers rij with n\ > ■ ■ ■ > nk > 0) give rise to the Nielsen integrals: @k,p{x) — Li fc+11P -i (x) (the 
subscript P -1 stands for p — 1 times 1). Thus, the representation can be rewritten as 

j-i 

Li,(l -x)=J2 ( Li j-r(l) - l.i.u •:,'•:■) (-If Liir(x), 

r=0 

wherein the term with Lij^l) should be dropped. Proving the asymptotics of the third sum for p = 2 
entails two curious representations of Nielsen integrals, l|l(J.9|l and H1U.10|) . which themselves boil down 
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to the following bizarre identities for Bernoulli numbers: For n a positive integer, n 2, 

(-!)>-£ " (-ly/n-A _ 1 
^ r ^ I \n-l) n 2 

r V f-^ \ « / r n - rJ n 

r—l l—r v 7 

Proving the asymptotics for higher p, one gets a further such identity, and a whole tower can be built up. 
The first identity is easy to prove, but the second has resisted our best efforts (and those of experts). 



2 Formal sums of Bernoulli numbers and zeta-values 

We start with formal infinite sums involving Bernoulli numbers. The notation is formal because the 
sums diverge [Bik ~ (2k) 2k ). Nevertheless, they are useful as one can recognize constant terms from 
the expressions ^ffyr c/c in the Euler-Maclaurin formula. We shall use such (diverging) expressions to 
recognize constants in future applications of the Euler-Maclaurin formula. 



Lemma 2.1. In formal notation: 



«2) = f + J> 3l (2.3 

«3) = 1 + ]>> 21 2 '^ (2.4 

<(4) = f + £^ (a + ?. ( r + 2) (2.5 

CW= 2^2+^(2l)! (i-1)! ( 6 

\ log(27r) = -C'(0) = 1-^2 J2k^t^l) < 2 " 7 
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(2fc)(2fc- l)(2fc-2) 
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LfJ 

-Bg+i rr A// _\ 1 \ " 5 2 fc gj (gg-g q -2fc + i) / i\qr | \ " f^2fe ,~ . . 

-j+T*q - (, K-D - -fi+Ty? - 2^ jwy. ( q -2k+iy. K' 1 ) 9 ! 2^ T^j (2k- a- 1) 1 

fe=i fe^n i+i " ' 



E#=7-| (2-12) 



( 2k 



E oT~\ = C(0) + 7 + \ = ~\ log(27r) + 7 + § (2.13) 



2fc- 1 
i 



E = -2C'(-1) + 2C'(0) + 7 + % (2.14) 



2fc-2 

2 



E = 3C ' ( ~ 2) ~ 6C ' ( - 1} + 3C ' (0) + 7 + 1 ( 2 - 15 ) 



2fc-3 

2 



E ^ = E(- 1 ) i+ ^fiV^- J ') +7+ ^-5-Eii J 1 ^) ; (?>o) (2.i6) 



! = 1 



Ssfe) = S H, 1'-i) ( ' (i - i) + T ^ (2 - 17) 



LIJ+i 



E B »— 2— =<(3)-<(2) + i (2.18) 

1 m B 1 

fl„=log„ + T +--X:^- 5E + 0( 17i l TT ) (2.20) 

<^ 2: > - «') - <d>FT + ^ - % m " l.-V " + °<^> < 2 - 21 ' 

Proof. The first five lines are the constant terms in the asymptotic expansion of H n , , , ff^ 4 - 1 , Hn ■ 
The next five lines are the constant terms in the asymptotic expansion of E k q log fc (q = 0, 1, 2, 3), the 
first being given by the Stirling formula (for logn!). These are the generalized Glaisher constants [F-03], 
see lemma 15.21 The third set of lines is obtained recursively by partial fraction decomposition from the 

previous: ^- = {2k} J 2k _ ]} - E£o ^L-P i while we used (2fc)(2fc-j) = "jCst ~ 2^7) for CEO- 
Lines (|2.18|l and i|2.19|l are miscellaneous linear combinations that we shall use. The last two lines yield 
the asymptotics for the generalized harmonic numbers Efc=i 'F' w ith C(*) m lEOH a s it appears in the 
Euler-Maclaurin formula. □ 

Lemma 2.22. 

El ^ (M+iza = C(i ) _ i±i (2 . 2 3) 
E ^ {2k -%^ l - 3) = C(0 - C(< - 1) + (i _ ^ _ 2) (2-24) 

E (2fc - j) -^y- 2 - j) = E(~ir (;)C« - r) + ( . _ _ 1} , (Kj<i-2) (2.25) 

= E(- 1 ) r (^C(i-r)+c,- i , ^i-l,i>2) (2.26) 
= Ef" J " 1 + rN )C( i - r )+ c ^ Cj<-l,i>2) (2.27) 



r=0 
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with C(l) standing for Euler's 7 in the last two equations. As to the constants, they are computed recur- 
sively, \/j 6 Z,i 2: 

c i,i = C j-M — Cj—i^—i, hence Cjj = (— ^ 'Y^r^i—^-Y Cj-i-\+ r , r ~ Cj-i+2,2^, which is useful for the 
fourth line; while Cjj = Cj\fi »_i — Cj,j_i, hence Cj t i = JX=3 c i+i!f — c j, 2 > * s useful for the fifth line. The 
third line can be used for values of Cjj with 1 ^ j ^ i — 2. For any j : Cj t 2 = ^5- ■ 
For instance: Cj_2,j = - ^ — (i 3, and c 2 — ~\), cmd c^u = (—l) l ^ 1 H i _i (i 2). 



Proof. These are direct combinatorial consequences of 1)2.6(1 . □ 

In the lhs of the fourth line (j > i — 1), one can choose to exclude the few non-zero terms at 
low values of fc; in that case '^2 1 is replaced by X)ri=i] an d the only change lies in the constants: 
Cj.2 = 2 Sfe=i ^2fe(l — and the recursions remain unchanged. 



3 Stirling numbers and their asymptotics 



For future use, we set 



S r ,s,t 



E 

s^ii<"-<i r <t 



'1 



which is for t < r + s — 1. Define So, s ,t :=lift^s — 1 and S r ,s,t := if r < 0. These numbers relate 
to the Stirling numbers of the first kind [ 1 , defined by ^2 k [^]x k := x(x + 1) • • • (x + n — 1), or to the 



signed Stirling numbers s(t, r), defined by ^ fc s(n, k) 



x 



x(x — 1) • • • (a; — n + 1), in the following way: 



S r ,i.t = [*+i] A- = (— l) r+ *s(t + 1, r + l)/t\. The generating function for these three versions are: 



I<lo g( l +.))•• = I>,^ 
r! 1 — ' LrJ n! 

1 T n 

±(-]Dg(l-x)) r = Y i Sr-l,l, n -l — . 



Here is the relation between Stirling numbers of the first kind and harmonic numbers. Denote by 
{r} = r^ 1 . . . r\ l a partition of the integer r into I different parts, ie r — i\r\ + • • • + i;rj. Then: 



^,i,» = (-i) r EIT 



ij\ 



(3.1) 



M j=i 

The first few cases are S\ t \,n — H n (partition {1}=1), and: 

S2,x, n = ~k H n ] + l H n partitions {2} = 2, l 2 . 

S 3 ,i, n = |^i 3) - \H^H n + \Hl partitions {3} = 3, 21, l 3 . 

Sa,i,u = ~\H^ + \H^H n + \{H^f - \H^H 2 n + partitions {4} = 4, 31, 2 2 , 21 2 , l 4 
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This can be used, in combination with (|2.20() and l|2.21[l . to compute the asymptotic growth of the 
Stirling numbers to arbitrary order. The terms contributing the most are those with highest power of 
H n (/~ logn): partitions {r} = l r ,2V , etc. Thus the asymptotic expansion starts as 

S r ,i,n = h (log n) r + ^ (log nf- 1 + (log n)^ 2 + . . . for fixed r (3.2) 

For an alternative proof of this result and for shedding light on the decreasing sequence of logarithms, 

see end of sectional Equation jjQt is the main result of [W-93] and was strengthened in [H-95]. 

As for the asymptotic behaviour when r grows as quickly as n, say for n — r fixed, formula 1(3. 1|) 

is helpless; but we can easily find the solution by intuition: n\S ni i tn — 1, nlSVi-i.i.n = 1 + ■ ■ ■ + n = 

n(n+l)/2, n!5„_ 2 ,i,n = E£=i E£=i 1+ i *i*2 = \ (square - diag) = ±[n 2 (n+l) 2 /4-n(n+l)(2n+l)/6] = 

n(n + l)(3n 2 -n- 2)/24 ~ n 4 /8. Similarly: n\S n - 3i i tn = E" 1= i E" 2 =i 1+ i E" 3 =i 2 +i *i*2«3 = |f (cube - 

plane) ~ h{n{n + l)/2) 3 ; and in general we will have n\S n -k,i,n ~ n( n ( n + l)/2) fc , that is: 

1 1 /n 2 \™~ r 
S'r l n = — r 7 \7 ( — ] + • • • for n — r const. 



n! (n 

These two asymptotic growths agree with the results of [MW-58] obtained by saddle-point evaluation of 
the generating function integral (a method already used by Laplace two centuries ago for Stirling numbers 
of second kind). The same results were re-obtained in [KK-91] from recursion equations using the ray 
method from optics. Formula 1(3. 1() . however, gives as many terms as desired for the growth with r fixed. 
Formula 1)3.1(1 can be inverted to yield 

4 Asymptotics of sums involving (logk) p /k q 

Lemma 4.1. ELi ^ - (logn) [H n \ logn - 7 ] + 7i + EZU + °(^)' 

Proof. Write the lhs as ^f^- + Efc=i ^T^' ^ 0T f( x ) : ~ ^aP> nere are ^ ne ingredients we need: J f(x) = 
/(2*-i)(a;) = (2k - iy H ' 2k -^ losx . Thus the Euler-Maclaurin formula tells us that J^lZl i2 r M = 
~ \ ^ + E (^3fc-i - log n) - E • The log n terms yield - (log n) E jagfrr , which 

estimates (log n) [if„ — 7—^7] by (|2.2U|I . Writing (log n) [£/"„ — ~ log n — 7] assures us that the remaining 
terms are inverse powers of n (easily tractable under the asympfe trick). The constant 71 is by definition 
the first Stieltjes constant. □ 

Using the same method of proof, we easily generalise. 

Lemma 4.2. £Li ^ = ^ + 72 + |^ " E^i %T ^- 2H ^^ +H ^~ H ^ + 
0( ( j^S), with 72 =E . 5tmtfori», for p > 0: 

V ( lQ g fc ) P _ 1 f| nn7 tf+i 1 -,. 1 i ( lo S n ) P y ELo^(log") P " r , n r (i°gn)^ 

fc - J+l^S 71 ) + >+2 ^2fe +U{ S n * m +l J, 

fc=l ' fc=l 



with d p . k , r = ^ \p- r )\ s rj,2k-i and -f p = (-l) p pi 2~2 ^ Sp,:L,2ft-i. 

fc>p/2 

Here, 7 P is the p-th Stieltjes constant by definition. For p = 0we have So. s ,t = 1 and 70 = 7 (though 
one needs to add ^ to j p in this case). Note that in the last sum, we could drop the requirement k > p/2, 
as S Pt i : 2k-i vanishes for k = 1, . . . , p/2. We will neglect this in future. 

Note that the expressions we find for j p are exactly the same as those one finds when directly an- 
alytically expanding the zeta function via the Euler-Maclaurin formula. One uses f(x) :— x~ s with 
/(2fe-D( x ) = - s ( s + 1) . . . ( s + 2k - 2) a;— a*+i = -(2k - 1)! Ep'o 1 S P ,i,2k-i(s - 1)? x-*- 2 ^ 1 : 

N M 

as) = lim E " = 7 + 5 + E T^Tf s ( s + !) ' ' ' ( s + 2A; - 2 ) + crror ( s ' M ) 



fc=i 



1 r+s:(i:^--.)(«- 1 )' 



p^O fc>p/2 



We also easily generalise in another direction: 

Lemma 4.3. £Li ^ = -C'(2) + (bgn)[ff^ - C(2)] - ± + Efc=i + 0( w), 

—("'(2) = 1 — 2 B2k{H2k — 1)- Similarly, for q ^ 2; 

fe=l fc=l 

wift -C'(g) = ^ - E ^r^^p^a^ - ^-0- 

The meta-generalisation regroups the two previous results: 
Lemma 4.4. For q 2 and p 0: 

(lQgfc) P _ / iw(p)/^ (g-V+i( l0 g") P ~' , i (lQg») P ^ ELo rf P,9,fc,r(lQgn) P ~ r ^ (logn)" , 

fe=l fc=l 

„,,«, W , _ B2fc (2fc+g-2)! (-l) r p! r, 

WtUL Ctp,g,fe,r — ( 2 fe)! (g-l)I (p-r)! c 'r,Q,2S!+Q-2) 

and (_1)P C (P) (0 ) = + (-l)VEl ^^^^,2^+^-2 - 

Proof. Write the lhs as — h Efc=i ^°k^ ■ For f(x) :— ^ lo ^ , here are the ingredients we need: 

f ff T \ _ 1 V^P P ! (l°g^) P ~'' f(i)|V| - 1 (g+'-l)! yP (-l) r +'p! r, Ana- t~\J>->" fNR- 

) J\ x ) — ~ ( q -iy+i x <i-i l^r=0 fr=r)\ ' J \ X ) - ^TT)1 Z^r=0 (p-r)! ^r^.g+i-lW^J l iNtS - 

S 1 ... = for r > i). Now simply apply the Euler-Maclaurin formula. □ 

Application to numerics of the ( function. Note that this formula, together with the asympfc trick, 
allows a very rapid numerical computation of (q) (for positive integer q) , much more efficient than 
current mathematical softwares. Lemmas 15 . 31 or T6 . 21 provide a formula for negative q. 
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5 Asymptotics of sums involving k q (logk) p 

Lemma 5.1. For p ^ 1 we have: 

n p m B 2k V^P (-l) r pl q i\nan\P- r 

E(^) P ^E^(^«) P - r +Klog^+con S t-E 2fc(2fc - 1)L ^ 1 ^^^^ +0(^g 

k=l r=0 k=l 

with const = (-I)" C (p) (0) = (-l)fp! ( - 1 + Ei 2 fc(ffc-i) Vi,i,2fc-a) . 
Lemma 5.2. For g we Ziawe: 

n 9+1 rfi 

^2 k q log fc = ^— [(logn) - ^ + §n« logn + £ J^^L [(logn) + (ff, - H q _ 2k+1) ] n^ 1 

k=l q k=l 

m 

+ COnst +(-l)V Yl llk^k-q-l) + Q(^=?) 

fc=rn+i 

vmn const - c, ( gj + ?+1 fl, - (g+1)2 2^fe=i (2fe)l ( ? -2fc+i)! I ^ ?-Z^fc>[§]+i (2fe)...(2l- 9 -i) 
feeing the generalized Glaisher constant of \2.11)) . For odd q, it is understood that the last term of the 
sum /jfc— i • ■ ■ (with k = |~|] and independent of n), which equals ^ H q , should not be counted as it is 
already counted in const. 

Again, the meta-generalisation regroups the two previous results: 

Lemma 5.3. For q and p ^ 1 we have: 



5>«(b g *) p = « ,+1 E^w^ + 5« 9 (i°g«) p + E (E^^iog^r 

fe=l r=0 k=l r=l 

+ COnst + V ^=l^P,g,fc,r(lQgTl) P ~ r , ^(log-V 
/ ' ^k — q— 1 



ana a P:9: fc, r ■— (2fc)--.(2fc-g-i) (p-r)! 2^j=ov ^ "-V.i.g^r-j-i.i^fc-g^ 



„„Yft ^, , — B 2k ql pi o 

WtUI, ^p,q,k,r ■— (2fc)! (g-2fc+l)l (p-r)! °r, ? -2fc+2,g 

._ B 2k (-\Y+"+ 1 p\q\ ^q 
■p,q,k,r ■— (2fc)...(2fc-g-l) (p-r)! Z^j= 

and const = (-l) p C (p) (-<?) + f$rP ! <Sp,i,9 = { ( q +ly$ ~ £fc=i c P,«,fc,P ~ Eft^m+i d P,<?,fc,P- 

r £ i 

For odd g, if is understood that the last term of the sum EfcHi • ■ ■ (with k — |~|] and r = p), which is 
constant and equals ^ p\S p> i t q, should not be counted as it is already counted in const. 



Proof. Write the lhs as n q (logn)P + ElCi fc 9 (logfc)P. For f(x) := x <! (logx)P, here are the ingredients we 

(g+l) r+1 ^r=0 (p- 



need: / f(x) = j^L, E,=o and 



5t E^=o T^F)^ogxf- r q\S r , q - i+l , q for i < g, (NB : 5... = for r > i) 



^=^Y? r = ( " i r; ) r > ' (log^- r 9lE^o(-l) J 'gj,l, g gr-j-l, 1 ,i- g -l for * > 9, 



Now simply apply the Euler-Maclaurin formula. □ 
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6 Asymptotics of sums involving (logk) p /(n — k) q 

Lemma 6.1. Yl=\ ^ = i^gnf + 7 flogn) - C(2) + El-Li - C'(l -*))£ + O(^). 

Proo/. Write the lhs as log(n - 1) + YZ=i i^r- Thc integral of /(x) := is: - log(a;) log(l - 

f) - Li a (|), thus 1 /rVW = (logn) 2 - C(2) + 2Li a (I). Note also that /«(*) = il T ^ T r + 
(J-'P-CV ^o that [f^\ X )Yr = (2fc-l)![log(n-l) + E^ 1 

In total: 

(n - 1)3 2k V j 2k- j J 

Now use Ij2.17ll as well as the expansions log(n— 1) = logn— + ^ + . . . ) and ( n \y — J2i^o d-i ) W^r- 
There are nice cancellations so that only — C'(l — k) survives at power -\. □ 

Lemma 6.2. For p ) 1 « have: 

n-l , s p p 

E = (logn)^ 1 + 7 (logn)f + Ec p , r (lognr i - 

fe=l r=l 

- E E<w( w - - *)) jjg + ^t 11 ^) 

fc=l r=l 

wt/j c Pi r := fe^Fjf C( r + 1) an d dp.r.k '■— fa-rM ^ — i,l,fc-i- Hence, the constant term is (— l) p p!£(p + 1). 
Proof. Write the lhs as log p (n - 1) + y™I 2 With /(a;) := we have: 

/ /(x) = -log(l-£)(lo g!E ) P + £ L T^W^r+l(mOgx) P - r 
J r=l 

f(x) = log(n) log*(n - 1) + £ ^i^Li r+1 (l - A) log p - r (n - 1) - (-l) p p!Li P+ i(A) 

1 r=l 

P 

= log(n) log p (n - 1) - I ( ^ 1 (ei,,(i) - C(r + 1)) (logn) 1 " - (-l) p p!Li P +i(^) 

(i°g^) p , „y» (-i)Vj ^ ( -ir,i 



(B)] »-i = (2fc _ 1)!log P (Ti _ !) + (2fe _ !), ± log-"(n - l) 5 - 1 ' 1 ^ 1 - 

j' = l ^ ' r=l 7 J ' 

where we used the notation &i :T (x) := Yln>i ~ J ~~^r 2 ~^ x n from section^] In total: 

£ (!2i^ = r _1 / ( x)+ iiog p («-i) + ( 7 - j)io g "(n-i) 
fc=i ™ fc 7i 

(n - 1) ^ 2fe V ^ <f- r > ! v ' 1 v ; 2k- 1 J 

±... + Jzl>l V ^i(y l r ^lo g P-'>-l) S '- 1 ^- 1 -(-!)> +.., 
(n - 1)J ^ 2fc V 4> v ' j K ' 2k- j J 



1 We have used: Li2(l — x) + Li2(x) = — log(s) log(l — x) + f(2) (proof by derivation). 
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The remainder of the proof are nice cancellations, which are impossible to prove in the general case; we 
exhibit here the case p = 2 as a pattern for all other cases. For p = 2 we have: 



E = r 1 + * log2 (" - +E ^[/ (2fe - 1) wir 1 

™ — K ./ 1 <JK 



(logn) 3 + 2(logn)[©i,i - C(2)] - 2[6i, 2 (±) - C(3) + Li 3 (i)] + i log 2 (n - 1) + 

2 v^Bsfe/, , , H 2k _ 2 \ 2 „ B 2fe / log(n - 1) Hi , ff 2fc -3 ^ 



W ^ (n - 1) Y 2fc V 2fc - 1 / (n - l) 2 ^ 2fc V 2 2 2k -2 J 

+ 2 (- 1 ) J V ^ 2fc / l°s(" ~ jj ^fc-j-i N 

+ (n - l)o fs 2k \ J J 2fc — j _/'**' 

L^J+i 

Now replace E L i J+ i % by 7 - ± - ErM and set h 3 := 2£ L|J+1 ^^fef 1 - The ^ ~ D wiU 

(-1)' 
»>I i(n-l) 4 



cancel out due to log(l — ^) = — log(l — = Eijsi Jn-i)' ( anc ^ tnc squared version of it). Now use 
(110. 9|) and i|l(J.l(Jfl to simplify the rest and arrive at 

E = (logn)3 + 7(logn)2 " 2 <( 2 >( 1o s") + 2 ^ + E [^r^ 2 ^" - ^-1) - Jr - E (-i) r (* Z J) h r] ^. 

Use C"(-9) = ~~ ( q +iy-i _ "i+r 25 ' 2 > 1 .9~I]iii ^fffr ( 9 _ 2 fc+i)! ^S 2 , q -2k+2,q -J2 k ^ rf 1+1 w^^k-l-i) (H 2 k- q -2-H q ) 

1 _ v-vi O-i)! 1 

2fc(2/c-j) ~~ Z^r=l (j-r)! (2fe) ■ ■ ■ (2/c-r) 



from lemma as well as the partial fraction decomposition 2 k(2k- 1 = EJ=i F-rll (2fc)---(2fc-r) ^° sn ow 



that 

» e t -»-»'(;:!)!<-&-')+*] 

or equivalently: — Er=i(~l) r (r-i)' 11 ' = — ^) + p'- This completes the proof for p — 2. The proofs 
for p ^ 3 run similarly. □ 

Lemma 6.3. For q ^ 2 we have: 

n—X 1 7 g— 2 m 



i=l fe=l 



with C q := 7 — iFj + ^. Hence, there is no constant term. 

Proof. Write the lhs as log(n - 1) + E^? With : = we have: 

/ /(s) = _^_ aogiB) ^ I W^_lo 8 (n-x)_g 

./ J W 9-1 'i'- 1 fa - n)9-! /g-1 ^ 



g — 1 \ a;'- 1 (ic-n)?- 1 / g-1 n?- 1 ^ i{q — 1) na-*- 1 (n — x)* 

1 t^ 2 1 / 1 \ 2 log(n - 1) 



■ log(> - 1) + 



g-1 ^ (ij-l)i n'- 1 -' V(n-l)' / g-1 n' 



1 1 ' tn-Ui ( n - r.\<t+i + 2s 



[/— ' wir 1 = m„ - 1) + (2fc - 1), e ( 2fc - r + r 2 ) ( — ^ 

(o — 1) r V o — l / V (n — 



( 1 


n q-l-i \ 


(n 


- I)' 






,-i ; 


(n - x 


9+ 




2k — 1 


(" 


T\ 




r V 



„i r V 9 _ 1 'M«-l) r (n - l) 2fc - r +'J- 1 ■ 

In total: 

logfe _ /•»-! i , , S 2fc (g+2fe-2)! 



E = + J lo s(" - 1) + log(" - 1) E 



+ 



(2fe)l (g-1)! 

_ V ^/ 2fe + 9- 3N l _ 1 1 y Ihk /2fe + g-4N 1 (-1)" Safe /2fe- 1\ 

1 ^ 2fe V g-1 / (n - l) 2 2 2-, 2 k V g-1 / "' (n - l)?- 1 g-1^ 2feVg-l/ 

1 



(n-1) 
(-1) J + 1 



(-!)"+ 1 B 2h ,2k -2s y, B 
2s 2k V g - 1 / V 



(n - 



g z ^2fcVg-l/ ^ (2fc)(2fe - 1) 
(-1)" ^ Safe /2fc- j -2\ p' + g- In „ 

q + j 2s 2k { q -l ) \ q -x ) 2s ( 2 fc)(2fe-j-l) 



in 



Now use <|2.17[1 . <|2.25[1 . as well as the expansions log(n — 1) = \ogn — + + . ..) and rf^zjxd — 
Ei>o ( J i"-7 )n?FT- There are nice cancellations so that only — C'(l — &) survives at power □ 

Note that without going through the proof, one can empirically determine the values of the C q just 
using the asympfc trick: one first uses the trick to quickly determine the 30 first values of C qi then uses 
it again to determine the asymptotic growth of those values up to 0(-4-) and recognizes the growth of 
harmonic numbers. 

Again, the meta-generalisation regroups the two previous results: 

Lemma 6.4. For q 2 and p ^ 1 we have: 



k=l ^ ' i=l r=l r=0 



(logfc)P ^ 

- E (fe©) E ^Oois«r - (-D p (^T 2 )c (rt (i - *) ) 

fe=l ^ r=l ' 

With Cp_i^ r '.— ^(p— r )! ^r— 1,1, i— 1 

_ (-r-fp! J - (P - r)S r ,l,q-2 + S r -l,l,q-2 7 + 21=2^-8,1,9-2 C( s ) for r = 0, 1 

ana «p,(j,r : — (p— r)! x ] 

[ EI=r-«+2 S'r-a.i.g-aCCs) for r = q,...,p (in case p ^ q), 
wherein F) rq are the rational numbers 

p — 2q — p + r c , „ q-j-2 „ r „l, * \i C C \ C 

V ir.lj-1 / S P -2-r,l,,-i-2 ^ ^ <=0 ( ~ 1 J 6i . 1 ■ r- 1 b P - 3 - r- j , 1 , 9 - j - 2 - r \ Op-l,l,g-2 „ 

D r,( := 2^ 2_- : : 2_- 7 — r~\ M ^ 6 p> 1,9-2- 

In particular, there is no constant term in the asymptotic expansion. 
Proof. Write the lhs as log p (n - 1) + Y^Zi with / 0) ■= we have: 

p 9 — 2 



■/ (9-1)! -fe, 3 'ni" 1 



^(p-fc)! I ^ \ ^ (i - r - 1)!(? - j + r)!i n«-'(n - ^'-'^ ^ (, - ljn*" 1 fe3+u »' 

p! / * (-1)' • 



P-2 ! /-.Nr + p 



+ 



E ( E *-^-V„i.«- 9 io g -(„ ^ d) Lij>-i(i - i))] - [(-i) v E (^i^-i ^-^ 1 ^)] 



(() = (g-l + i)! (logx)P * p!(logx)P-'- (-i)jj!( g + i-j + r )! 

(9-1)! (n f^, (p - r)\(q - 1)! ( j + , 



(g — 1)! (n-a;)9+ i ^ (p — r)\(q - 1)! ^ (j + r ~ - j - r + 1)\ (n - x)i+ i - r -i + 1 x r +i- 1 

,.(»*-!) M1 »-l _ (g + 2fc-2)!logP(»- 1) (2fc-l)! 2 ^ 1 1 ^ (-l)^- r p!logP--(n-l) ( g + 2fc-j-2)! 

1 (9-1)! (9-1)! ^(n-lK^ (p-r)! j(2fe-l-j)! r1 ^ 1 

(2fc-l)! 2 ^' 1 1 (-l)P+ipl(q + j -l)! o 

(9-1)! f^ Q (n-l)«+i (2*-j-l)il 

The total expression for EsUi (n-fcjg f rom the Euler-Maclaurin formula is too messy to write out. As 
usual, cancellations will be hard at work and the result will boil down to the rhs in the lemma. The 
closed expression for the D r q was particularly hard to find (empirically). □ 
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Application to the asymptotics of Stirling numbers. 

Had we not known the asymptotic growth of Stirling numbers (|3.2|) , we could easily find it by induction 
from the leading terms in lemma FOl and RT2l Assuming the empirical result (via the asymp^ trick) that 
the coefficient of x n in (— log(l — x)) p has leading behaviour ~ p(logn) p_1 /n, we prove: 



(p + 1)! £ S P ,i tn -i— = (- log(l - x)f+ 1 = (- log(l - - 
^ — ' n *■ — ' n 



?i * — ' rt 

n-l „ , n p _i 



g /y (lQgfc)^ 1 | y (logfc)* 



fe=i fe=i 

( + (log nf) = (p + 1) V 

n n 

hence <Sp,i, n ~ |f(logn) p . The same inductive proof works for the next-to-leading term of (|3.2|) . In that 
case, we also need the next-to-leading term of lemma IB. 21 For each subsequent term that we want to prove 
in (|3.2|) . we need one more term of lemma RT2l while the leading term of lemma B"l2l is enough. One thus 
sees how the sequence of decreasing logarithms in (|3.2() is intimately related to that of Er=i c p-r(log n) p ~ r 
in lemma W7R 



7 Asymptotics of sums involving l/k q {\.ogk) p 

Lemma 7.1. For p ^ 1 we have: 

n—l ^ p— 1 ^ ^ m 2k — 1 ^ ^ 

^ (togT^ = ^ h(?l) ~ n ^ (ioi^f + ~ 2 (bg^ ~ 5 ( 5 (logTj'+p) + 

with c PjT := an d d p , r ,k '■= 2fc(ffc-i) ^ p ^i)\ ' ^r-i,i,2k-2, an d Cp,o « s i^e constant term. 

The log-integral is defined by li(z) := J 2 

Proof. This follows from the Euler-Maclaurin formula with fix) :— (\ Q ^ x y and 



/ m = (An w<» - * E 



(r-l)l 



Si (P - !) ! (log a:)*' 



i - 1)! x 4 (loga;) r +!> ' 

From these, it is also straightforward to write down the 'exact' expression for the constant C Pj o> involving 
a formal (infinite) sum over Bernoulli numbers. We omit it as it is not enlightening. □ 

Note that the second sum on the rhs is just the start of the asymptotic expression of the first term, 

since li(n) « nJ2 r >i jisfffi • So we mi S nt replace the two terms by ^lyr E r > P (lo^ny ■ Trlis is indeed 
what one obtains when numerically looking for the asymptotics of the lhs; the first term is (± g n \ P , and 
correctly so. Yet since this asymptotic expansion diverges for all values n, the replacement would be 
disastrous for numerical evaluation of the constant C Pi q. 
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Lemma 7.2. For p ^ 1 we have: 

n—l 1 1 1 _. m 2k— 1 , 1 

E 1 _ g ^ 1 1 ( "p.^ fc ^ 1 l of i ) 

„A:(logfc)f pa (p- lHlogn)?- 1 2n(logn)P £j V £jj (logn^+W n 2fe l "" m+ij ' 



fc=2 



wif/i d p ^ r _ k := ■||-^^jp5'r,i,2fe-i7 C P ,i * s ^ e constant term. 

For p = 1, the second term on the rhs has to be replaced by log(logn) (which becomes the leading term). 
Proof. This follows from the Euler-Maclaurin formula with f(x) :— - ^ gx y an d 

/ /(a) = - ( P -l)(lo 8a) P-i 



/ W (z) = (-!)**!£ 



^ (p-1)! s' + i (loga;) r +f 



Lemma 7.3. Forp ^ 1 and q ^ 2 we have: 

p-i 



□ 



1 1 

fe ^fc9(logfc)P ' (P-i)! i/-o-v n9 _i ^ (l gn) r 2 n« (log n)P 



m 2fe~l 



S ( S fW n)r+P ) n 2fc-l+g + (^T?)' 



(log n) r +PJ n 2k - 1+ i 

fe=l r=0 v to ' 

withcp. q . r := (1 ~' ?)i | p _ 1 ) ! (r ~ 1)! anddp^k ■= ^i^-iy' ^{p-t)]' S r, q ,2k+ q -2, andC VA is the constant term. 
The exponential integral function is defined by the principle value of the integral: Ei(x) := — f™ x ^j-dt. 
Proof. This follows from the Euler-Maclaurin formula with f(x) := tt— ^ and 

./ (p-l)l a:*" 1 ^ (p-1)! (logx)'' 

/0 (x) = ( _ 1) i(*+lzi)if : (P-l + r). 



(9-1)! ^ (p-1)! xi+i (\ogxy+p 

□ 

Again, the third term on the rhs is just the start of the asymptotic expansion of the second term, since 
Ei((l-g)n) « (l-g^iogn)'- - So we mi g ht replace both terms by the infinite sum n^ 1 J2 r > P (log »V • 

But since this diverges for all n, the replacement is disastrous for numerically computing the constant 
c 

^p,q- 

Note that when p ^ 2, the previous lemma makes sense also for q = 1, and one recovers the preceding 
lemma (since (1 — q) p ~ 1 ~ r vanishes unless r — p — 1). 

For large p or large g, it is quite obvious that the main contribution to the sum X)fe=2 fcg (logfc)p comes 
from the term k = 2 and that the constants C p . q will converge towards 2q ( log2 ) P ■ Just how quick they 
converge can be empirically determined: asymptotically for large p or q, we have 

n 1 L p-ap-bq , -cp-dq , 

p ' 9 2«(log2)P + + + '--' 

with a = 0.09405,6 = 1.0986, c = 0.3266, d = 1.386. Of course, these values are nothing but log(log3), 
log 3, log(log4), log 4, so as to obtain g^j^gp + 49(^4^ ! So we come back from where we started. This 
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comes as no surprise when C p<q = X)fc^2 fcg (i og fc)p > ^ u t ^ * s a sur P r i se when the infinite sum does not 
converge, ie. when C Viq is not the leading term in the asymptotics, eg. when q = and p becomes large. 





5 — u 


i 
i 


9 


Q 



i 

-t 


n — 1 • 
p — 1. 


— n 94394 


n 7Q4fi7Q 




fl 937QQR 




2 : 


3 1 0399 


9 10Q74 


699606 


305808 


143463 




4 96079 


9 06589 


889388 


41 991 4 


1QQ0Q1 


4- 


00344 
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Zj . rjrjjy ± Zi 


1 1&Q9& 

1 . XOc'ZiO 








7 4fi^74 


3 49989 




8086^9 




fi- 


9 Q9ni 5 


4 75831 


9 33093 


1 1 51 06 


^71 944 

U . U 1 -L Z 44: 


* V lu 6 


q — u 




2 


3 


4 


p=l: 


1.4427 


0.721348 


0.360674 


0.180337 


0.0901684 


2: 


2.08137 


1.04068 


0.520342 


0.260171 


0.130086 


3: 


3.00278 


1.50139 


0.750695 


0.375348 


0.187674 


4: 


4.3321 


2.16605 


1.08302 


0.541512 


0.270756 


5: 


6.24989 


3.12495 


1.56247 


0.781237 


0.390618 


6: 


9.01669 


4.50835 


2.25417 


1.12709 


0.563543 



Table 1: Comparison between C P:q and 2 g(i 1 



We may want to add li(2) = 1.045 to Ci t o so as to obtain the constant 
lim„^oo(£!£2 - f 2 l j^) = 0.80192543. Similarly, we add log(log2) = -0.36651 to Ci,i, so as 
to obtain the constant lim n _>oQ ( £fe=2 fc(iog fe] ~ I2 a:(iogx) ) = 0.4281657. Both values already occurred 
in [B-77], see also [F-03]. We were not able to recognize an exact form for either of these two constants 
(using PARI for integer linear combinations of other constants, or using Plouffe's inverter or his Maple 
code). 



8 Asymptotics of the Taylor coefficients of (— z/log(l — z)) k 

We now use lemma 17721 to generalise a result known to Polya [P-54] about the Taylor coefficients of a 
certain generating function. In 1954, Polya [P-54] noted that 

1 z 

a » n — ^ for /W = ^^ — r : E fl » zn t 8 - 1 ) 

n(logn) z — log(l — z) 

We shall be interested in the asymptotics of the a n when the generating function is raised to some 
power k (positive integer). For k — 1, the series begins as 1 — \x — j^x — ... and all coefficients are 
negative except ao- The a n for k = 2 are asymptotically given by the convolution of those at k = 1, viz. 
a n = £ j(iogi)^ (K-i^iogtw-i)) 2 • Since this sum makes only sense for i running from 2 to n — 2, we write 
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the terms - n(logn) y ~ ( n - m $ n -i))* - ~ n(logV twice separately. By symmetry, we can write: 

1 (n-2)/2 

2 



n(logn) 2 ^—^ i(logi) 2 (n — i)(log(n — i)) 2 ' 

where in the last sum, log(n — i) logn/2 and = + From lemma IV. 21 we know that 

i(logi) 2 = C ~ (logn/2) + ^( (logn/2) 2 

_ r rru-.i re: i.„ tx._j.i 

(n — i) (log i) 2 

(log l) z v (log 7 

2C 2 



X] n ' 2 t(lo g t)2 = C — (i og »/ 2 ) + Q( (io g n/2) 2 )' wherein the constant C is figurative, since the quantities 
l/(i(logi) 2 ) only approximate the exact values of the Taylor coefficients. Further, ^"/ 2 (n _ i) j losi)2 ^ 
^E n/2 7lSgl F -0( 7I5 ^ F ). Overall: 



n(logn) 2 rt(logn) 2 rt(logn) 3 Vn(logn) 4 / 

Since the singularity of f(z) 2 at z = 1 is of higher order than that of f(z), the decrease of coefficients 
should be stronger; hence the n ( logn p terms have to cancel each other and so C = 1/2. We are left with 
a ™ ~ ~ n(iogw) a + • • • • One similarly obtains: 

>( 1 rj-rz ) for f( z ) = (— — Z - r) . (8.2) 

\n\oen) k+2 J w V-Wl-zW v ; 



n(logn) fe+1 \n(\ogn) k+2 J V — log(l — z) 

A naive attempt at justifying Polya's result (|8.1|) would be to use Cauchy's formula a n = ^ § c ^t+i z 
and to compute the contour integral on the unit circle, z = e l8 . Note that — log(l— e l8 ) = — log(— 2ie l6 / 2 sin 
Thus we would have (wrongly) 



I f 27r e -i(»-l)fl rf g ^ 1 yc/n ^ ^ c/(2tt) / 1 

27rio - log(2sinf) ~ 2vr Jo -log0 ~ nlogn Vn(l 



f 1 ) 

Vn(logn) 2 / 



where we replaced J v by J Q C since for large n, only small values of will contribute substantially to 
the integral. This approximation, however, does not yield the desired result - presumably because one 
cannot replace e - l ( n - x ) e by 1. Similarly, had we used partial integration with f'(0) = 



we would have ended up with Jo L lo A® ea c ^^ n , again with the wrong leading term. As the 



(i^-log(2sinf)) 2 ' 

2TTin Jo 0(log0) 2 1-0 n log n ' 

integral is not tractable by the Laplace method, the saddle point method or any other trick described 
in [dB-58], we shall see in the next subsection that the solution lies in a clever choice of the contour of 
integration. 

The result (|8.2(l is not new, but was already obtained by Norlund in 1961 using combinatorics of 
Bernoulli polynomials, and rederived by Flajolet and Odlyzko in 1990 by evaluating the contour integral 
in Cauchy's formula. For completeness, we present these two alternative and elegant paths below. 

8.1 The Flajolet-Odlyzko approach 

In 1990, Flajolet and Odlyzko summerised the 'transfer properties' of analytic functions, viz. the be- 
haviour of the function at the first singularity on the convergence radius is directly reflected in the 
behaviour of the Taylor coefficients. One of their result is [FO-90]: 
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Theorem 8.3. (Flajolet-Odlyzko, 1990) Let f(z) be analytic in \z\ < 1 + rj except for a singularity at 
z — 1, and let 

/(z) = o((l-z) Q (-log(l-z)r) as z->l (a,7GR). 
Then the coefficients a n in f(z) — ^2 a nZ n grow like a n — O ( "j ^ . 

Proof, (sketchy) . It is comforting to see that the proof boils down to a mere application of the Cauchy 
formula, ie. a contour integral around the origin, viz. a n — § c ^*} + f z , but the contour has to be 
chosen cleverly - as in figure ^ 




Figure 1: The contour of integration, excluding the singularity at z = 1. 

The contour C will be a circle of radius 1 + rj with a tiny roundabout around the singularity; the main 
contribution will come from this little near-circle C\ (of radius 1/n) around z = 1. Note that for any 
compact domain inside our contour C, there is a constant if such that \f(z)\ ^ K |(1 — z) a (— log(l — z)) 7 |. 
On the circle C\, we have 1 — z = e z9 /n and the following bounds: |/| ^ K(^) a sup| log(ne _l6 ')| 7 — 
^ l °^ ^, as well as |z|™ +1 ^ (1 — — ► 1 and J Ci \dz\ ^ 2ixjn. Hence the main contribution to 
the contour integral can be estimated by: a n » ^ § Ci ^^l+f^ = ( ^a"l J , as we wished. □ 



This result was readily obtained, but is treacherous when a is a non-negative integer, say 0: the 
Taylor coefficients of — log(l — z) decrease like 1/n and not like (log n)/n, and those of l/(— log(l — z)) 
decrease like l/(n(logn) 2 ) and not like l/(n(logn)). For this case, it is useful to have a precise asymptotic 
development, whose derivation we sketch as follows (see [FO-90] for details). Let f(z) be the function 
(1— z) a (— ~ log(l— z)) J . Change to the variable z — 1+t/n and expand log(— rt/i) 7 = (logn) 7 ^i>o (J)( — 
io ^~^ y ■ The contour integral now contains the piece G.j := ^ J C2 (— ^) a (l°g(— t))' l e^ t dt, whose contour 
can be deformed to a well-known integral that simply yields Gi = d l a Y(- a ) • Hence the main contribution 
from the Cauchy formula for f(z) — ^2 cinZ n is: 

^-^PE 1 ^^ f - /(z) = (1-^-^(1-^. (8.4) 

The first term of the sum, i = 0, yields r(-«) an< ^ thus simply drops out when a is a non-negative integer 
(where r(— a) has a pole). This explains the above treachery. For a — and 7 = —k (negative integer), 
we have ("*) = (-l)^*** -1 ) and g t := ^ a | or ^y = -1, 2 7 , ~3 7 2 + (i = 1,2, ...), and so obtain 



Hi 



our sought-for asymptotics: 

{ i )9i 



n (logn) I+1 n(logn) 2 



for 



(8.5) 



8.2 Norlund's approach 



Most surprisingly, Ij8.4|l was arrived at 20 years earlier from quite a different angle, namely by Norlund [N- 
61] who recognized the Taylor coefficients a n of f{z) as being values of generalised Bernoulli polynomials: 



(1 - *)«(-- log(l - z)V = ]T [ -J-Bi n+ ^(a). 



(8.6) 



n>0 



The Bernoulli polynomials of order 7 are defined by 



and coincide with the usual Bernoulli polynomials for a — (and thus yield the Bernoulli numbers for 
7 = 1). To arrive at Q8.6JI . use Cauchy's formula ~ B "„ ! < ' a ' > = § t ^ e t_ 1 y r rft, then substitute i = log(l — z) 
and shift 7 — ► 7 + n + 1; you thus obtain the Cauchy formula equivalent to (|8.6|l . At a = 1 and 7^0, 
we recover the Stirling numbers of first kind: 



(-1)" B^ +1 \l) = 
The polynomials satisfy 



n + 7 
7 



n + 7 

7 



7117! 
n + 7 



"7—1,1,71 



+7-1 



B { r?\a)= B% +V >{t)dt and (a) = (a - l)(a - 2) • • • (a - n) 



(-!)> - a) ..- (1 - a) = (_ ir r (n a + 1) 

r(i - a) 



Hence we can readily obtain an asymptotic expression for a — and 7 = — 1: 
(-l)"B< B) (a) f 1 Y{n+\-a-t) 



r(l-a-*)r(n + l) 



(It 



J vr(i-« 



— a — f ) 



1 e -tlogn 



9y 



-rft + Of-r— ) 

1 Vn log n / 



.(logn)r(l -a) J -logn T(l-a-t) 

where we used Stirling's approximation n\ ~ e nl ° STl ™ + 2 lo s( 27m ) m ^ ne se cond step and partial integration 
in the third step. Further partial integrations yield 
(-1)" B^ia) n a 



where the derivative term evaluated at a = 1 is: 9' 



y 1 

fWlogn) i + 1 a T(l- 



alir(l-a) L/ alor(-a) 



(!"«) 
91 



0( 



fori = 0,1,... . Overall: 

(-1)" si n) (i) 



E 



^ rt(logn) I+1 



for f{z) 



(logn)''+ 1 ) 

: 5 i,with fli = 0,-l,2 7 ,-3 7 2 + 

z 

-log(l - z)' 



17 



From here, it is straightforward (via induction) to deduce the general case 7 = — k (negative integer) 



n\ j^ i n(\ognY+ 1 n(logn) 2 JK> \-\og[l- z)J 



in full agreement with (|8.5|1 



9 Asymptotics of the incomplete Gamma function 

This sections examines whether the constant C\ % \ met in section 7 also occurs in the constant term of the 
the function Y^n°=2 ( n \o\n) B • answer is Y es : as the general theory shows, but 7 has to be subtracted 
to obtain the full constant term. The resulting expansion is presented in (|9.2(l . and on the way we shall 
derive the following intermediate result about the asymptotics of the incomplete gamma function (a is 
an arbitrary constant): 

POO 

T(-s,as):= / t~ s ^e~ l dt = -log(as) - 7 + s[±(log(as)) 2 + a - ci] + 0(s 2 log 3 (as)) as s -> 0. 

J as 

9.1 Preliminaries 

Recall the coincidence in the constant terms of the following asymptotic expansions: 

fc^i s 

n-l 1 1 

Z) «l°gn + 7 - — + ... (n^co) 
fe=l 

Landau confirmed this coincidence for a broader class of Dirichlet series: suppose ^2 n<x h(n) ~ ax + . . . 
(among other constraints on h(n)), then: 

ft.(n) a 



f * rj s q — I 

for some constant /?, and 



n s s — 1 

n— 1 



rs alogx + p + . . . (x — > 00). 

n 



n 



We shall be concerned with a weaker generalisation. First recall the discrete partial integration formula 
for some continuous function <fi and some sequence a n with primitive A(x) :— Y^ n <x °™ : 

y a n (f>(n) = A(x)(f)(x) — / A(x)(f> (x)dx. 

n=a+l Ja 

When a, b G Z and a n = 1 with A(x) = [xj , the formula reduces to 

A b fb 

> 4>(n) = x4>(x) — / hxJ^'Oc)^ + ^( a ) 
a •/» 

= [<!>+[ (x- Lxjy + 0(o). 
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If (j> s (x) is a suitable function depending on a parameter s, like 4> s (x) = the derivative in the second 
integrand will ensure that we may exchange the limits b — > oo and s — > 1 (since a; — [^J is bounded). In 
other words, the first integral contains the singularity as s — > 1, while the second integral yields merely 
a constant. Denote by (f>\ the function <^> s obtained after taking the limit s — *■ 1; we then perform the 
partial integration backwards: 

/b b-l „k+l 6-1 /.fc+l b r-b 

(x - [zj)# = £ / (x - fc)0i = E U(fc + 1) - / 0i) = E - / 

/c=a J ' £ k=a Jk k=a+l J a 

Having previously taken the limit b — > oo, we obtain our desired generalisation: 

b ,b 



lim (J2Mn)- <M = Hm ( $>i (k) - / 
For instance, for s (a;) = and a = 1, we have: 

lim fc(s) — = lim fj? 6 - logb) = 7. 

s— >1 V S — 1/ fc^oo V / 

For </>(x) = ^jog^-p and a = 2, we have by lemma fT"2l 

OO -. /"OO 7 ^ 1 

lin l ( E n — - / ( 1 X \. ) = „ lim ( E h lo g( lo s 6 ) + lo s( lo s 2 )) = c ^ + lo s( lo s 2 )- 

s->i V z — ^ (nlogn) s J 2 [xlogx) s / b^oo V z — ^ nlogn / 

(9.1) 



9.2 Gamma function asymptotics 

In order to find the constant term in the asymptotics of the function 'Y^ =2 ( K i gn) s » we s ^ u need to 
expand the corresponding integral up to the constant term. This will include finding its singularity at 
s = l. Note first that 

f°° dr f 00 
/ ?xl^ = (S - 1) r(1 - S ' (S - 1)1 ° g2) = (S - 1) / , t ~ Se ~ tdt - 

J 2 {XLOgX) J(s-l)log2 

The incomplete Gamma function is defined by T(z, s) := J s °° t z ~ 1 e~ t dt. Hence = r^^| + T s . To 
simplify matters, we shall first study the behaviour of T(— s, s) as s — > 0. Note that T(0, s) ~ — logs — 
7 + . . . and r(— s, 0) ~ — - — 7 + . . . , so that the simultaneous limit will behave like the weaker of both, 
ie lim s ^or(— s, s) ~ — logs + . . . . Equipped with this intuition, we proceed by expanding d s T(— s, s) 
(around s = 0) and then integrating the expanded result. Now with z(s) = — s we have 



-fr(-s,s) 

as 



dT dz dT 
dz ds ds 



r s - x e _t (logt)dt. 



The first part is expanded as —e~ s s~ s ~ 1 = — Ie~ ;5 ( 1 + logs ) = — 1 + (logs + 1) + . . . and gives us - upon 
integration - the required (— log s) term, so that the second part can contain at most log-singularities. 
Differentiating the latter gives us e~ s s~ s ~ 1 log s + t~ s ~ 1 e~ t (log t) 2 dt, which behaves as (log s)/s + . . . , 
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so that the original integral behaves as ^(\ogs) 2 + const + .... In general, J s °° t s x e t {\ogt) l dt will 
behave as — ^j-(logs)'+ const + . . . , and by bootstrapping we obtain: 



r^e^dt = -(logs) + c + s[±(logs) 2 + 1 - ci] + . . . 
t-' -1 e-*(logt)(ft = -±(logs) 2 + ci + s[±(logs) 3 + (logs) - 1 - c 2 ] + 



i- s - 1 e - t (log^di = -^ I (log S 



i-i+l 



^(logs) l+2 + ^ ^ii^ ( logs) r - c l+1 



+ 



In this way, we can start with the n-th line at order 0{s) and recursively determine the expansion of the 
first line up to 0{s n ) in terms of the constants of integration Co, . . . , c„. The latter can be empirically 
determined: c = -7, c x = 0.98905, c 2 = -1.81497, c 3 = 5.89038, c 4 = -23.568, etc. We have computed 
the first 300 of them and it seems that their asymptotics are Cj = (— l) i+1 (i! — e 8 llogJ + . . . ). It would 
be interesting to know more about these constants. In particular, we obtain our desired expansion: 



3 j_ £2 

2 



T(-s, s) = - (log s) - 7 + s[| (log s) 2 + 1 - Cl ] + s 2 [ - i (log s) 3 - (log s) + I + 
Next, study T(— s, as) for some constant a. Now j^T = Y z {— 1) + T s a. Bootstrapping yields now: 



(*-Q) 



t-^e-^logtydt = -^(logs)^ 1 + Cl + s[-^(logs) 4 + 2 + 



(-i) r+i »i 



(logs) r - c l+1 



r=0 



+ 



with the constants Cj taking the same values as before. In particular, 

T(-s,as) = -log(as)-7+s[i(log(as)) 2 +a-ci] + s 2 -±(log(as)) 3 -alog(as) + f--^+a +... (s -> 0) 
Thus we can answer the question above: 

f l + (s - 1) log(s -!) + ...][- log((s - 1) log 2) - 7 + (s - l)[i log 2 ((s - 1) log 2) + log 2 - Cl ] + . . 



(x log a;) 



= - log(s - 1) - (7 + log log 2) + ( s - 1) [ - I log 2 (s - 1) + (log 2 - cj - 7) log(s - 1)+ 

+ log2-ci] + 0((s-l) 2 log 3 (s-l)) 

Recalling H9.1|l . we deduce the constant term in the asymptotic development of the original function: 



00 ^ 

E n — ^ ~ " lo s( s " V + c ^ ~ t + °(( s - V lo e 2 ( s - x )) ( s "> V 

(n\ogn) s 



(9.2) 



n=2 



where Ci,i = 0.794679... was given in table □ 



10 A representation of polylogs and of Nielsen integrals 

In this section we present a representation of polylogs in terms of Nielsen integrals which we have come 
across while embarking onto the proofs in section 6. Conversely, in l|10.9|) and l|10.10|) we give repre- 
sentations of Nielsen integrals involving Bernoulli numbers, which boil down to mysterious identities 
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for Bernoulli numbers and harmonic numbers fconiecture 110. lT|) . We only treat the cases 6i t i(x) and 
61,2(1), but are convinced that similar formulae hold for all Si iP (x). 

Lemma 10.1. 



Li 2 (l - a;) = C(2) + (logx) E - - E ^2 ( 10 - 2 ) 



n * — ' n 

1 



™71 TT , XJ"(2) TT 

Li 3 (i -x) = c(3) - c(2) E - - a°g*) E % i * B + E (— + ( 10 - 3 ) 

Z ^ 77 Z ^ 77 Z ^ V 77. 7-7Z / 



n z — ' n * — ' v n 
1 11 

n-l Hi_ 



n tt yvn— ± flj-i 

m(i- a; ) = c(4)-c(3)E— + c(2)E^ zi ^ + ( 1 °g ;E )E 1 ' " 

Z ^ 77 Z ^ 77. Z — ^ 



X — 



n *■ — ' n * — ' n 

1 1 

t=l V i 1 ~ 1 , 



_ y- / L, i= i i + p + L,i=i * \ x n ( 10 4 ) 

' — ' V n n / 

1 

Li, (1 - z) = C(i) - C(i - 1) E + C (j - 2) E ^f^V - • ■ • + (-1)^(2) E 



1 



(-^•(log^E ^" 2 ' 1 '"" 1 ^" - E f ^" 2 -"" 1 + 5j ' 24 2 "' 1 )x" (10.5) 



n * — ' \ n n 
1 1 



with Tj_2,n-i '■= ^tttjt 3 fc2 (sum oner a// 1 ^ i\ < ■ ■ ■ < ij-3 ^ n — 1 and l^fc^n — 1, fc 7^ 
ii, . . . satisfying the following recursion: Tj_2, n -i = ^17=1 Tj ~ + "^"^'a''^ 1 ■ 77ie Zasf siim 0/ 

^10.5p can also be written as ~^2iJ^&i,j-r~i(x) l ° e |.|~" C ' ) wif/i 6i iP (x) := Ylm^i ^ £ ~^r ] —^-x n , that is 2 
xd x &i. p (x) — - — — for p 1. T/ims the above can be rewritten as (j ^ 2) 

0,(1 -a;)=E (C(i - r) - ©w-r-xW) (10.6) 
wherein the term with £(1) should be dropped. 

Proof. Show the recursion A Lij(l — x) = Lij_i(l — x) using the fact that 

i:iZ l (9 x &i, j -r-i(x)) } ^^l = WZ^£l j-, 1 )!; 1 ;;, = (binomial formula for (1 - l)^ 1 ). Note 

that log could be replaced by Lzj£iig2) or by (— l) r Liir(x), where the generalised polylog is defined 

by 

Ex" 1 
— 5 j--. Thus the formula could read: 
n-i . . . n t fc 

ni>--->n fc >0 1 * 



Lij-(1 - a) = E ( L V-r(l) - Li 2jlJ -.- 2 (x))(-l)'- Li lr ( 

r=0 



wherein the term with Lii(l) should be dropped. 



□ 



The Si jP are special cases of so-called Nielsen integrals (for k 1): 



z We have used (— log(l — a;)) p = p!S n >0 [p] nT = P'En>0 Sp ~ 1 ^ L '"~ 1 which is proved by expanding: lhs 
dr\ (1 — x)~ r = . . . (notation and convention from [GKP-89]) 
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so that &k-i,i(x) = Lifc(cc) and xd x &k,p{x) = *5k-i,p(x)- Also: ©1,1(2;) = 1^2(2;), ©1,0(2:) = log x. 

As a corollary, equate the last term of p().4(l with the corresponding quantity in 1)10. 6[1 and find: 
n E^l* + T^) = E^i 1 (ipNj* ~ ! or, using partial fractions: 

n-l ri-1 tt(2) „ 

(n-i) 2 V j i 2 

»=i v ; i=i 

which can easily be generalised: 
Lemma 10.7. For p 1: 

n-l „ n-l P-l rr(p-r+l) „ 

(n - i)P ~~ ^ V ^ i r 

t=l v ' i=l r=l 

Proof, lhs = E^i* = X^Ji 1 EjJr 1 IFJ = E"=i Eti 1 i?cFIJ = rhs > usin S the P artial fractions 

1 _ v p 1 1 1 n 

Lemma 10.8. Wii/i 6j := Er=i ^ awe /o r < 1- 

eW.j-UM.Wfc^-^)') (10.9) 

<^(*> - E ^ - E ( ^g^ - + Li l( ,)](-^) 1 ) (10.10) 

n^l i>l ^ 7 

Proof. The first equation is equivalent to an identity for Bernoulli numbers: 



r * — ' / \n — I J n 

r—l l—r 

in which the second sum equals I Y^WG) = Ei"=o (-!)'(?) = -^C-i) = "^O- 
Thus the equation boils down to the well-known identity for Bernoulli numbers: E^=o (™)^ r = ^ n or 
mnemotechnically: (B + 1)™ = £>„ (upon replacing B k by P&). 
The second equation is equivalent to 

r ^ I \p - I J 1 1 2-^i \ 2-^i r 2—, 1 \ n - I 1 1 p - n p 2 

r=l l=r V ^ 7 n=l v r=l i=r V / / ^ ^ 

or 

^ (-1)^,. * (-1)' /p- 1\ ^ ^ (-1)»P„ | | 2 H p _! = g p _i 

r I \P ~ v ■ / -^?i 2 (p — n) p p 2 p 2 

Note also that 



where the second equality follows by induction, while the first comes from 

ms = (-if- 1 y (-if -h^hi = (-1Y- 1 la - X y +1 ~ log(1 ~ ^ 1 =(-ir- i y ( ~ 1)P 'Y p 

1=0 1=0 v ^ 

Hence the second equality is equivalent to the following conjecture, which has resisted the author's best 



efforts. □ 
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Conjecture 10.11. For p a positive integer: 




P 



1 



r 



)=H?l 1 + 1 -H, 



p-l- 



Conclusion 

Though an application of the Euler-Maclaurin formula is nothing distinguished, it turns out that its 
use for the four sums of sections 4,5,6,7 resp. brings along a wealth of by-products about Stirling 
numbers, their relation to harmonic numbers, their asymptotics, about mathematical constants and 
their representation as formal (diverging) sums over rational (Bernoulli) numbers, about more general 
asymptotics of complex functions (incl. incomplete gamma functions), as well as algebraic manipulations 
on polylogs and Nielsen integrals. All this research was only possible because we used the asympfc trick 
and numerical mathematics. 

Extensions of this paper would be doing the same for other sums containing logarithms, foremost 
Y] fc 9 /(logfc) p , ^2(n — k) q (logk) p , but we do not expect any new property. Perhaps only more unknown 
mathematical constants would come to light. Otherwise, the constants C Pj9 of section 7 still await an 
exact form. Further, one could attempt to prove conjecture IIP. 1 ll or write down the similar (and more 
complex) identities that one obtains when carefully going through the proof of lemma l6~2l Perhaps even 
more bizarre identities would show up by carefully analyzing what happens in the proofs of lemmas 
and 
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A The asympfc trick 

Assume we are given numerically the first hundreds of terms of a converging sequence (s n ), (n £ N), and 
that its asymptotic expansion goes in inverse powers of n, ie. s := cq + — + % + . . . . 
Goal: determine the coefficient cq. 

Trick (by Don Zagier): apply the operator ^\d k n k on s (k G N) to find 



Hence this gives k more digits of precision for Co, as long as k is not too big (ie. the binomials not too 
big). Call this operation asympfc. In practice, the operator d is the difference operator As :— s n+ i — s n . 



co + (-l) 



k °k+l 
n k+l 




+ 
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To determine c\, subtract Co, multiply by n and apply asymp^ or: differentiate (ie. take successive 
differences) and multiply by — n 2 . 

The crucial point in the success of this trick is that the errors generated by the difference operator 
on a monomial are themselves monomials of lower powers: A~ = f^ppp — -~ = — ^ + J ^"j+^ 2 + • • • . 
These will be swept away at the next applications of A. The same would not be true if the operator 
acted on terms like (log n,y . 

NB: we can see whether the k decimals are correct by checking the convergence of the series: if we 
have 400 terms, say, write every 80th term in a column (ie. 5 terms in total) and see how quickly the 
digits agree. 
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